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Abstract 

We consider a chain of nonlinear oscillators with long-range interaction of the type l/l 1+a , 
where I is a distance between oscillators and < a < 2. In the continues limit the system's 
dynamics is described by the Ginzburg-Landau equation with complex coefficients. Such a system 
has a new parameter a that is responsible for the complexity of the medium and that strongly 
influences possible regimes of the dynamics. We study different spatial-temporal patterns of the 
dynamics depending on a and show transitions from synchronization of the motion to broad- 
spectrum oscillations and to chaos. 
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A chain of nonlinear interacting oscillators is a model of the wide-spread in- 
vestigation of different physical phenomena such as synchronized behavior of the 
system, bifurcations to different regimes, spatial-temporal turbulent or chaotic 
dynamics, different instabilities, appearance of defects, and many others. The 
long range interaction between oscillators leads to a new qualitative dynamics 
and thermodynamics. Our consideration is related to the power-law long-range 
interaction that makes it possible to find new regimes of the system and to es- 
tablish a link of the corresponding equations of motion to the equations with 
fractional derivatives. 



I. INTRODUCTION 



The goal of the paper is to consider different dynamical regimes, from synchronization 
to chaos and turbulence, in a chain of large number of coupled nonlinear oscillators with 
long-range interaction (LRI) of a power type. The potential of interaction is nonlocal and 
proportional to l/l 1+a with I as a distance between oscillators and < a < 2, (a 7^ 1). It 
will be called a-interaction. Transitions between different regimes of the chain behavior are 
considered as a function of a. In the continuous limit the system reduces to the fractional 
generalization of the Ginzburg-Landau (FGL) equation and the chain of oscillators can be 
considered as a discretized model of FGL, called DFGL. Since we also consider the model 
with complex coefficients, it either can be related to a fractional generalization of the complex 
nonlinear Schrodinger equation. 

The literature related to this type of problem is fairly wast and we would refer only to 
some reviews and closely related articles. The complex Ginzburg-Landau (CGL) equation is 
typically considered for pattern formation in different media IlL \A yfl. CGL equation appears 

nnnd 

in numerous physical models (see for example [4|, |5j, |6j, 17]). Long-range interaction with finite 
radius of interaction was considered for complex media in j^, and for the a-interaction in 

;he chain of coupled oscillators 



jiol . 11. 12. I3I . Different regimes of synchronization i n t 



161 ] while the synchronization 



with nearest neighbour interaction can be found in reviews [15 

with a-interaction for DFGL was found in {l?! fisj ] . The FGL equation was introduced in [l9( 
to describe the wave propagation in complex media, when the dispersion law has fractional 
power of the wave number. The related case appears in description of weak turbulence . 
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More examples related to the chains with LRI and their solution and breather type solutions 



and applications can be found in 



21 



22]. 



The main part of this paper is a numerical simulation of the chain of TV = 512 coupled 
oscillators with a-interaction and constant external force that pump energy into the system. 
This system links to the FGL equation in continuous limit based on the formal procedure 
introduced in [l7, 23]. 



It will be shown that the forcing oscillatory media with LRI exhibits transition to chaos 
and turbulence when a decreases. As another interesting regime, we found a "collectivized" 
limit cycle that has a broad spectrum in phase space of any individual oscillator. Such a 
cycle exists within an interval of a. 

In Sec. 2, we present the basic equation for coupled oscillators and their continuous 
media limit. This part includes some modified version of the results of 



17 



24). In Sec. 3, 



we consider fractional generalization of the CGL equation and some of its stability criteria. 
The Secs4-6 play an auxiliary role to identify an interesting domain of parameters to be 
studied. The main numerical results are presented in Sees. 7-9. It is important to note 
that the main studied object is the chain of nonlinear oscillators, and a part of the obtained 
results can be appropriate for the complex FGL equation at least for a finite time. The 
FGL equation can also be used for some estimates of the critical parameters to change the 
regimes of behavior of the chain of oscillators. 



II. CHAIN OF OSCILLATORS WITH LONG-RANGE INTERACTION 



Following |24j, let us introduce the system by its action 



S[Z n 



-Z n (t)Z n (t)-V(Z n (t)) 



J2 U(Z n (t),Z m (t))), (1) 



n,m=—oo 



where Z n are displacements of the oscillators from their equilibrium, and U(Z n (t), Z m (t)) is 
defined by 

U(Z n (t), Z m (t)) = ^g J a (\n - m\) (Z n (t) - Z m {t)f, (n ^ m) (2) 

with go as an interaction constant, and with interparticle interaction 

1 



J a (\n - m\) 



\n — m\ a+1 



[a > 0). 



(3) 
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The usual linear oscillator term Z%(t) can be included into V(Z n (t)). The corresponding 
Euler-Lagrange equations are 



dZ n 
~dt 



(4) 



m=— oo 



where F(u) = dV{u)/du. 



Continuous 

to Z(x,t) 
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imit of equation (j4]) can be defined by a transform operation from Z n (t) 



23]. Firstly, define Z n (t) as Fourier coefficients of some function Z(k,t), 



ke[-K /2,K /2},i.e. 



Z(t,k) = Z n (t)e- lkx "=F A {Z n (t)}, 



(5) 



where x n = nAx, and Ax = 2n/K is a distance between nearest particles in the chain, and 

I f+K /2 



Z n {t) 



dk Z(t, k) e ikXn = Tl l {Z(t, k)}. (6) 

J-K /2 

Secondly, in the limit Ax — > (K — > oo) replace Z n {t) = (2ti/ K )Z(x n ,t) — > Z(x,t)dx, 
and x n = nAx = 2irn/K — > x. In this limit, Eqs. (Q, Qj are transformed into the integrals 



Z(t,k) 



dx e~ lltx Z(t,x) = T{Z{t,x)} = lim F A {Z n (t)}, 

Ax^O 



i r+oo 

Z(t, x) = — dk e ikx Z(t, k) = ^{ZOt, k)} = lim F?{Z(t, k)}. 
2tx Az-o 

Applying ([5|) to (j4]) and performing the limit ([7]), we obtain 

dZ(t,x) d a Z(t,x) 



Of 



+ 9o 



d\x\ 



F(Z(t,x)) = 0, (l<a<2) 



(7) 
(8) 

(9) 



where d a /d\x\ a is the fractional Riesz derivative 25| defined by 

d a Z(t,x) 



d\x\ a 



f- l {\k\ a Z{t,k)}, 



and 



/ixa\ 

g a = 2g (Ax) a T(- a ) cos (-) 



(10) 



is the renormalized interaction constant. Equations of type with different nonlinear 
terms were considered in 
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231 ] . For other values of a one can obtain similar equation to 



j9]) by performing the corresponding transform operation. 



III. COMPLEX NONSTATIONARY GINZBURG-LANDAU EQUATION 



Consider a chain of the nearest-neighbour coupled nonlinear oscillators described by the 
equations 

-jj_ z n(t) = (1 + ia)Z n - (1 + ib)\Z n \ 2 z n + (ci + ic 2 )(Z n+1 - 2Z n + Z„_i), (11) 

where we assume that all oscillators have the same internal parameters. Such equations can 
be obtained from (Tj0) for a = oo and the corresponding choice of go and V(Z n ). Transition 
to the continuous medium assumes that the difference Z n+ \ — Z n is of the same order as Ax, 
and the interaction constants C\ and c 2 are fairly large. Setting c\ = g(Ax)~ 2 , c 2 = gc(Ax)~ 2 
with new constants g and c, we get 

—Z = (1 + ia)Z - (1 + ^)|Z| 2 Z + #(1 + ic) — Z, (12) 

which is a time-dependent complex Ginzburg-Landau (CGL) equation. 

Let us come back to the equation for nonlinear oscillators with fractional long-range 
coupling 

^-Z n = (1 +ia)Z n - (l + ib)\Z n \ 2 Z n + g 1 (Z n - Z m ), (13) 

at *—f \n — m\ a+1 

where Z n = Z n (t) is the position of the n-th oscillator in the complex plane, 1 < a < 2, 
and the corresponding choice of F(Z n ) = (1 + ia)Z n — (1 + ib)\Z 2 \Z n has neen done. The 
corresponding equation in the continuous limit is 

-Z = (l + ia)Z-{l + ib)\Z\ 2 Z + g{l + ic)—Z, (1< a < 2, a ± 1) (14) 

where g(l + ic) = g a a , 

a a = 2r(— a) cos(7tq;/2). (15) 

Eq. (j!4p is a fractional generalization of the complex nonstationary Ginzburg-Landau (FGL) 
equation f[T2"j) . Equation (|T5|) is the same as (0J with a special choice of V. It was suggested 
in [3] to describe complex media with fractional dispersion law. 



The FGL equation (THl) can be presented as the system of real equations 
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Jt Y ^('-'W^t'- 9C 9^> ~ + Y2)(Y + bX) ' <16) 
where X = X(t, x) and Y = Y(t, x) are real and imaginary parts of Z(t, x). 

The results of numerical simulation of the system ( fl3l) can be applied for a finite time 

for the FGL equation (114j) although it is well known that typically the discrete systems are 

more chaotic than their continuous counterpart. The system (fl3l) will be called discretized 

FGL equation or simply DFGL. 

IV. PARAMETERS OF STABILITY OF PLANE- WAVE SOLUTION 

In this section, we obtain some conditions of instability that will be implemented in 
numerical simulations of the system (fl3j) . These conditions are easy to derive and interpret 
considering the FGL equation. 

To obtain a particular solution of the FGL equation ( fT4l) with a fixed wave number K, 
we consider Z in the form 

Z{x,t) = A{K,t)e iKx . (17) 

Substitution of (EE?} into (P3J gives 

^-A(K,t) = (l + ia)A- (1 + tb)\A\ 2 A - g(l + tc)\K\ a A. (18) 
The plane-wave solution of (|T8|) is 

A(x, t) = (1 - g\K\ a ) l l 2 e iKx - iUa{K) \ 1 - g\K\ a > 0, (19) 

where 

u a (K) = (b-a) + (c-b)g\K\ a , 1 - g\K\ a > 0. (20) 
Solution of (TTTJT) can be presented as 

X Q (x, t) = (l- g\K\ a ) 1 ' 2 cos (Kx - u a (K)t + 9 ) , 

Y (x,t) = (l-g\K\ a ) 1/2 sm(Kx-uj a (K)t + 6o), 1 - g\K\ a > 0, (21) 

where X = X(x,t) = R,eZ(x,t), Y = Y(x,t) = lmZ(x,t), and #o is an arbitrary constant 
phase. These solution can be interpreted as synchronized state of the oscillatory medium. 
In fact, it will be shown by simulation that the synchronized solution exists also for the 
DFGL equation jTJ]). 
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To obtain the stability condition, consider the variation of (fl8l) near solution (I2T1) : 

j t 5X = A 1X 8X + A 12 5Y, j8Y = A 21 8X + A 22 8Y, (22) 

where SX and SY are small variations of X and Y, and 

A u = 1 - - 2X (X - Wo) - (X 2 + F 2 ), 

A 12 = -a + gc\K\ a - 2Y (X - bY ) + b(X 2 + Y 2 ), 
A 21 = a- gc\K\ a - 2X (Y + bX ) - b(X 2 + Y 2 ), 
A 22 = l- g\K\ a - 2Y (Y + bX ) - (X 2 + Y 2 ). (23) 
The conditions of asymptotic stability for ff22l) are 

An + A 22 < 0, A U A 22 - A 12 A 21 < 0. (24) 

Substitution of Eqs. ([ZID and ([23]) into gives 

A n + A 22 = -2(1 -g\K\ a ), (25) 

A n A 22 - A 12 A 21 = (b(l - g\K\ a ) -(a- gc\K\ a )^ (36(1 - g\K\ a ) -(a- gc\K\ a )). (26) 



Then the conditions (|24j) have the form 

l-g\K\ a >0, 

1 - g\K\ a < a/b - (c/b)g\K\ a < 3(1 - ^^1"), (27) 
i.e. the plane- wave solution ffl9l) is unstable if the parameters a, b, c and g do not satisfy 



Condition (I27p defines the region of parameters for plane waves where the synchronization 
can exist. 

In the numerical simulation, we use the parameters 

a = -1.2, b = -2, c = 2, g = 1. (28) 
Then Eq. fT27|) gives the inequalities 

< 1 - \K\ a < 0.6 + \K\ a < 3(1 - \K\ a ). (29) 
As a result, the plane- wave solution with E = is stable for 

0.2 < \K\ a < 0.6, (1 < a < 2). (30) 



In our numerical simulation, we use the initial conditions with 

\K\ = 2tt/64 w 0.09817 (31) 

that is K is in the unstable region for 1 < a < 2. In the following the initial conditions 
with perturbation E ^ will be out of the boundaries (I3"0l . and evolution of the initially 
unstable states will be studied. 



V. MAPPING THE NUMERICAL DATA 

Numerical results are obtained as solutions of the coupled equations 
d 1 + ' N 1 

-Z n = (1 + ia)Z n - (1 + l b)\Z n \*Z n + £ (^ - Z m ) + £, (32) 

where n = 1, .., N, Ax = 1, and E is a constant external force. In simulations, the number 
of oscillators was N = 512, and the external force was E = 0.3. In the continuous limit, Eq. 
( 1321) transforms into the FGL equation with forcing. 

For all sets of parameters, we integrate the DFGL equations with the initial conditions 

Z n (0) = A e lKn , (33) 

where \K\ = 2n/T, the space period is T = 64, and A is an initial amplitude. Numerical 
solutions were stored at each t q = 0.005g, where q G N. 



To visualize the numerical solution we consider the following values. 

(1) We plot the surface \Z(x,t)\ 2 and the phase-space projection of the central oscillator 
(n = 0, x n = 0) formed by the variables: 

A(t) = \Z(0,t)f, Mt) = ^- (34) 

(2) In addition to ([34]), we plot Y(t) = Im(Z(0,t)) vs. X(t) = Re{Z{0,t)) for better 
resolution of the central oscillator behavior. 

(3) We calculate the discrete Fourier transform of the sequence Z(0,t q ) defined as 

Z n (uj) = — X] Z n (t q ) exp(-vi ujj tq), (35) 

^ q=Q 
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O-i 

Znitq) = ^ Zn(uJj) exp(i UJj t n ) , (36) 
j=0 

where ujj = 2irj/Q, (j = 0, ...,Q — 1). The power spectrum S of the sequence Z n {t q ) for 
q — 0, Q — 1 is given by 

^S(^-) = l^)| 2 - (37) 
Our main goal is to compare solutions of the DFGL equation for different values of a G (1, 2) 
and consider emergence of chaotic dynamics of the chain of oscillators as a function of a. 
The larger is a, the weaker is long-range interaction. 
In the simulations, we consider the following plots. 

(a) Color plots present surfaces \Z n (t)\ 2 vs t and n; 

(b) Plane (A, dA/dt) displays a projection of the trajectory of the central oscillator (see 
definition in (J34J)); 

(c) Plane {Re Z,Im Z) displays a projection of the complex amplitude Z = Z(0,t) of the 
central oscillator ClS db function of time; 

(d) Plots (log 

10 u -> l°gio S) describe the spectrum of time oscillations of Z(0, t) (see definition 

in Eq. flST}). 



VI. SOME NUMERICAL RESULTS FOR CGL EQUATION 



In this section, we provide some numerical results derived in 17|, |26j, |27|, |28|, |29| for the 
complex Ginzburg-Landau (CGL) equation in order to compare them, obtained for a = 2, 



with our results for a < 2. For many other results and details see 



15| 



(a) In paper [27|, the CGL equation has been considered for the parameters 

a = 2, a = 6= 1.333, c = -1. (38) 

and the phase turbulence has been observed. 

(b) It was shown in 28] that for the parameters 



a = 2, a = 0, b 



•1.4, c = 0.6 



(39) 



solution of the CGL equation has chaotic states and it is spatial-temporal intermittent, 
(c) It was found in 



281 ] that for the parameters 

a = 2, a = 0, 6=1.2, c 



-0.6 



(40) 



solution o: 
(d) In 



: the CGL equation has zigzagging holes near the transition to the plane waves. 
261 ] . the CGL equation was considered for the parameters 



a = 2, a = -1.2, b = -2, c = 2. (41) 

For E = 0.35 the solution displays pitchforks without defects while for E = 0.23 the rare 
defects were observed. 



(e) In 17|, the FGL equation was considered for the parameters 

1 < a < 2, a = 1, b = 0, c = 70. (42) 

It was shown that the solution had one stable fixed point for ocq < a < 2, (qq = 1.51) that 
corresponds to synchronization of oscillators. Decreasing of a below leads to a limit cycle 
via the Hopf bifurcation. 



VII. TRANSITION FROM SYNCHRONIZATION TO TURBULENCE 

Regular propagation in time of the initial state of the chain of oscillators will be called 
synchronization. To consider a-dependence of transition from synchronization to turbulence 
near a = 2, we use the parameters 

a = -1.2, b = -2, c = 2, A = 0.2, (43) 

similar to (14T1) that were used in 



26| for the CGL equation (fT4"]l . but for the chain of 
oscillators with a / 2. The simulation was performed for the chain of 512 oscillators with 
parameters equivalent to ( 1431) . 



In Fig.l., the surfaces of |Z n (t)| 2 vs t and x = n for a = 1.940, a = 1.935, a = 1.930, 
a = 1.925, a = 1.910, a = 1.810 are shown. For a = 1.940, we have a regular space 
structure of the plane-wave type. For a = 1.935, there is a space modulation, which deforms 
the regular space structure. For a = 1.930, we can see the appearance of defects and 
pitchforks. For a = 1.925, the structure in the space-time demonstrates sharp and drastic 
changes. There exist many defects and pitchforks. For a = 1.910, the number of defects and 
pitchforks is increased, and some pitchforks are joined. For a = 1.810, we can see chaos and 
turbulence, and synchronization is completely lost. We see that the amplitude turbulence 
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FIG. 1: Alpha-dependence of transition from synchronization to turbulence near a = 2. Surfaces 
of \Z n (t)\ 2 vs t and x = n for a = 1.940, a = 1.935, a = 1.930, a = 1.925, a = 1.910, a = 1.810. 
Simulations are realized for FGL equation with parameters a = —1.2, b = —2, c = 2, Aq = 0.2. 
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FIG. 2: Alpha-dependence of transition from synchronization to turbulence near a = 2. Plane 
(Re Z, Im Z) shows projection of the complex amplitude Z = Z(0, t) of the central oscillator as a 
function of time for a = 1.940, a = 1.930, a = 1.925, a = 1.810. 

is characterized by persistent creation and annihilation of pitchforks. The decreasing of 
order of fractional derivative means increasing a role of LRI. It is worthwhile to comment 
that the loss of synchronization and the emergence of amplitude turbulence is fairly sharp 
with fairly small change of a. 

In Fig. 2., the plane (Re Z, Im Z) shows projection of the complex amplitude Z = Z(0, t) 
of the central oscillator as a function of time for a = 1.940, a = 1.930, a = 1.925, a = 
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vlr — tre — rr " Ma 

FIG. 3: Alpha-dependence of transition from synchronization to turbulence near a = 2. Plane 
(A, dA/dt) shows projection of the trajectory of the central oscillator in phase space, where A(t) = 
\Z(0,t)\ 2 for a = 1.940, a = 1.930, a = 1.925, a = 1.810. 

1.810. For a = 1.94 and a = 1.93, there is a stable node, which means the existence of 
synchronization. Conformation of that can also be found in Fig.3, where 

the plane (A, dA/dt) shows projection of the trajectory of the central oscillator and 
A(t) = \Z(0,t)\ 2 . For a = 1.94 and a = 1.93, the attracting point that corresponds to 
synchronization of the chain of oscillators is clearly seen. Characteristics of the turbulent 
motion, shown in Fig.l. for a = 1.925, a = 1.91, a = 1.81, are different from the cases of 
larger a. This difference can be better recognized from the Figs. 2-4. Particularly, in Fig.2. 
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FIG. 4: Alpha-dependence of transition from synchronization to turbulence near a = 2. Spectrum 
of time oscillations of Z(0,t) (see definition in Eq. (|37p ) for a = 1.940, a = 1.930, a = 1.925, 
a = 1.810. 

for a — 1.925 and a = 1.81 behavior of ReZ = ReZ(0,t) and ImZ = ImZ(0,t) displays a 
disordered process that on the plane (A, dA/dt) in Fig.3. reveals a structure similar to what 
is usually observed for stochastic attractors. Nevertheless, more specific statement needs 
more detailed analysis since our system is open and has many degrees of freedom. 

In Fig. 4., we show the spectrum of time oscillations of Z(0,t) (see definition in Eq. (J37))) 
for a = 1.940, a = 1.930, a = 1.925, a = 1.810. For a = 1.94 and a = 1.93, the spectrum 
has small numbers of harmonics, and it corresponds to the regime of synchronization in 
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the oscillator medium. For a = 1.925, the spectrum is filled by additional harmonics. For 
a = 1.81, the spectrum has a dense set of frequencies for the region 0.1 < uj < 1 that reflects 
the chaos behavior and turbulence of the oscillator medium. For the region 1 < uo < 100, we 
can see the integer power-law for the spectrum S(u) ~ oj~ 2 that corresponds to the main 
frequency dependence of equation (ITB1 . 



VIII. TRAVELING WAVES AND BROADENING OF THE LIMIT CYCLE 



Large number of oscillators and periodic boundary conditions permit to consider localized 



traveling waves (see for example for discrete systems in 30|, |31()- Such waves were observed 



in our system (1521 near the values of a ~ 2 and disappear near a ~ 1. In the latter case, it 
will be also demonstrated space-temporary localized topological defects. 
The simulation was performed for the parameters 

a = 1, b = 0, c = 2, A = 0.2. (44) 

For the value 6 = 0, the nonlinear term is real. 

The waves propagation can be characterized by the group velocity v aj9 = dui a {K) / OK . 
From Eq. (I2"0l . we obtain 

v a , g = a{c-b)g\K\ a - 1 . (45) 

The phase velocity is 

v a , ph = oo a (K)/K = b —^ + (c - b)g\K\ a -\ (46) 

A 

For the parameters used in the simulations we have \v a , g \ > \v a ,g\- in our simulation 
K = 27r/64. For our simulation K = 2tt/64, and the decreasing of the order a leads to the 
increasing of the group and phase velocities. 



The corresponding simulation is shown in Figs. 5-8 for the surfaces |Z„(t)| 2 (Fig. 5), ReZ 
and ImZ plane (Fig.6), (dA/dt, A) - plane (Fig. 7) and spectrum S(u) in Fig. 8. Each of the 
the plate has the some additional information about solutions. There is a strong difference 
between values a = 1.99, a = 1.45, and a = 1.07, a = 1.05, a = 0.91 while the case of 
a = 1.09 is not clear since the finite time of simulation (t < 1000). The first case (a = 1.99, 
a = 1.45) shows traveling waves along x with regular (periodic or quasi-periodic) patterns. 
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FIG. 5: Broadening of limit cycle. Wave front deinclination. Surfaces of |Z n (i)| 2 vs t and x = n 
for a = 1.99, a = 1.45, a = 1.09, a = 1.07, a = 1.05, a = 0.91. Simulations are realized for FGL 
equation with parameters a = 1, b = 0, c = 2, Aq = 0.2. 

In Figs 5,6, we observe the approach of the solution to a limit cycle with a few harmonics 
in the spectrum (Fig.8). At the same time the smaller is a, the smaller scales (larger values 
of K) enter the solution. 

For a close to one (a = 1.07, a = 1.05, a = 0.91) a fairly irregular pattern of traveling 
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FIG. 6: Broadening of limit cycle. Plane (Re Z, Im Z) shows projection of the complex amplitude 
Z = Z(0, t) of the central oscillator as a function of time for a = 1.99, a = 1.45, a = 1.09, a = 1.07, 
a = 1.05, a = 0.91. 

wave declines at some time, and synchronized oscillations take place. In the phase diagrams 
in Figs 6,7 a broadened limit cycle type picture corresponds to collective oscillations of the 
chain with fairly rich spectrum presented in Fig.8. The closer is a to one, the shorter is 
time of break of traveling waves, and the synchronized breather type solution appears. Since 
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FIG. 7: Broadening of limit cycle. Appearence of topological defect. Plane (A,dA/dt) shows 
projection of the trajectory of the central oscillator in phase space, where A(t) = \Z(0,t)\ 2 for 
a = 1.99, a = 1.45, a = 1.09, a = 1.07, a = 1.05, a = 0.91. 

the growth of wave numbers K of the solution, its amplitude can reach zero giving rise the 
topological defects [3]. Their appearance is clearly seen from Figs 6,7 for a = 1.05 (see also 
Fig. 5), when the amplitude A(t) = |Z(0,t)| 2 reaches zero value. It is seen that for a = 1.07 
and a = 0.91 the dynamics is close to the appearance of the topological defect. Zero of the 
complex field Z result in singularity of the phase 9 = argZ. In two dimensions, points of 
singularities correspond to quantized vortices with topological charge 

n = ^-J> VM, (47) 

where L is a contour encircling a zero of Z. 
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FIG. 8: Broadening of limit cycle. Spectrum of time oscillations of Z(0,t) (see definition in Eq. 
Ij3711) for a = 1.99, a = 1.45, a = 1.09, a = 1.07, a = 1.05, a = 0.91. 



For 1 < a < 2, we have the limit cycles around the point A — 1, dA/dt = 0. It is easy to 
see the broadening of these limit cycles, when a decreases. 



Finally for this section, let us comment on the spectrum of time oscillations of Z(0,t) 
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for different a, presented in Fig. 8. The broadening of limit cycle leads to widening of the 
spectrum. The spectrum is filled out by different harmonics that are localized in the region 
0.1 < uj < 1. For uj > 1, we have the dependence S(uj) ~ u~ 2 that follows directly from 
equation (TH1) . 



IX. AMPLITUDE DEPENDENCE OF THE TRANSITION TO TURBULENCE 



In this section, we would like to show that the turbulent regime developing depends on 
the initial amplitude. For simulation, we use some parameters as in (l4"Tj) 261 ] except of 
a = 1.45. 

It is seen from Fig. 9, that the decreasing of the initial amplitude Aq from 1.3 to 0.005 
doesn't change, at least visually, the space-time structure of turbulence, but the change 
increases the time of developing an instability and beginning of turbulence. The turbulence 
possesses a robustness with respect to change of the amplitude of initial excitation and 
initial conditions. 



In Fig. 10., the plane (A, dA/dt) shows projection of a trajectory of the central oscillator. 
There are random rotations on the plane and a broadening of the spectrum in a finite region 
of frequencies. The Fig. 10 shows many loops of the trajectory for which dA/dt = and A 
is close to zero. That means that the turbulent regime includes many different points that 
are close to be the topological defects. 



X. CONCLUSION 



The main goal of this paper was to study influence of long-range interaction (LRI) on 
the developing of chaotic or turbulent motion in one-dimensional chain of large number of 
nonlinear oscillators. The LRI is characterized by the power of interaction a. In the con- 
tinuous limit the corresponding equation is the nonstationary generalized Ginzburg-Landau 
(FGL) equation with complex coefficients and fractional derivatives of order a along the 
coordinate variables. Our preliminary research show different interesting regimes of behav- 
ior of the chain depending on the value of a. We have observed a synchronized motion 
of the chain with different complexifications such as defects, chaos, space-time turbulence, 
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a = L45 f A - 1.3 a = 1.45, A = 

FIG. 9: Amplitude dependence of transition to turbulence. The decreasing of the initial amplitude 
from 1.3 to 0.005 gives that the space-time structure of turbulence almost is not changed. There 
is only increasing of time of begining of turbulence. 




a = 1 .45. Aft = 1,3 

FIG. 10: Amplitude dependence of transition to turbulence. Plane (A,dA/dt) shows projection of 
the trajectory of the central oscillator in phase space, where A{t) = |Z(0,t)| 2 , and the spectrum of 
time oscillations of Z(0,t) for a = 1.45 and A = 1.3. 

traveling waves. A possibility to use a continuous analog of the chain such as FGL equation, 
simplifies some estimates although it is well-known that discrete model is "more chaotic" 
than the continuous one. 
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